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Duality of Multiple Root Loci 

Hwangrae Lee and Bernd Sturmfels 


Abstract 

The multiple root loci among univariate polynomials of degree n are indexed by par¬ 
titions of n. We study these loci and their conormal varieties. The projectively dual 
varieties are joins of such loci where the partitions are hooks. Our emphasis lies on 
equations and parametrizations that are useful for Euclidean distance optimization. 
We compute the ED degrees for hooks. Among the dual hypersurfaces are those that 
demarcate the set of binary forms whose real rank equals the generic complex rank. 


1 Introduction 


Univariate polynomials of degree n correspond to points in a projective space P n . The 
multiple root locus A a associated with a partition A = (Ai,..., A d) is the subvariety of P" 
given by polynomials that have d distinct roots with multiplicities Ai,..., A^. The dimension 
of A a is d. The singular locus of A a is the union of certain codimension one subloci A M , as 
described in [[15], §3]. The degree of A a was determined by Hilbert in [Tl]. He showed that 


deg(A A ) 


d\ 

m,i\m 2 \ ■ ■ ■ rripl 


A1A2 • • • A d, 


( 1 ) 


where rrij denotes the number of parts A, in the partition A that are equal to the integer j. 

The multiple root loci A a have been studied in a wide range of contexts and under various 
different names: coincident root loci [5] [9] [15], pejorative manifolds [16] [23], strata of the 
discriminant mmm, A-Chow varieties [19], factorization manifolds [24] Definition 5.2.4], 
etc. Our motivation arose from the desire to understand the geometry of a model selection 
problem considered at the interface of symbolic computation [T3] and numerical analysis [23]: 
given a univariate polynomial h, identify a low-dimensional A a such that h is close to A a . 

Finding a point in A a that is closest to a given h is a problem of polynomial optimization 
[4]. We here characterize the geometric duality that underlies this optimization problem, in 
the sense of [3] Chapter 5]. The key player is the dual variety (A a ) v . This variety lives in 
the dual projective space P rt and it parametrizes all hyperplanes that are tangent to A a . 

The duals to multiple root loci were studied by Oeding in [T9]. He shows that (A a ) v is a 
hypersurface if and only if mi = 0, i.e. all parts of A satisfy A* > 2. In m Theorem 5.3] an 
explicit formula is given for the degree of the polynomial that cuts out this hypersurface: 


deg((X A ) v ) = A±hL.( Al -l)(A 2 -l)...(A,-l). (2) 

m 2 - * * ‘ m p- 
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For the application to optimization in [3} Theorem 5.23], this is the number of complex 
critical points one encounters when minimizing a general linear function over an affine chart 
of A\- For instance, consider n = 5 and A = (3,2). Following Example 14.21 and [B] §3], the 
polynomial for (Aa) v is the apple invariant of degree 12. So, optimizing a linear function 
over quintics with a triple root and a double root leads to solving an equation of degree 12. 

The present paper is a continuation of the studies by Hilbert El and Oeding |I9]. It 
is organized as follows. In Section 2 we set up notation and basics. In Theorem 12.51 we 
parametrize the conormal variety ConA that links Aa and (Aa) v . Theorem 12.91 offers a 
parametrization for the projective duals of multiple root loci. These results are derived from 
the apolarity theory for binary forms, as described in the book by Iarrobino and Kanev [12] • 

In Section 3 we study the multidegree of the conormal variety Coiia, and we summarize 
what is known about the ideals of Aa and (Aa) v . Table [0 offers a census of small instances 
up to n = 7. We also examine duality for hooks, and inclusions among the dual strata. 

Section 4 offers an application to tensor rank over M. Theorem 14.11 characterizes the 
algebraic boundary of the set of binary forms whose real rank equals the generic complex 
rank. When n is even then this involves the use of Chow forms [TO] and Hurwitz forms [2T] . 

In Section 5 we turn to Euclidean distance (ED) optimization. Theorem 15.11 determines 
the ED degree of Aa for hook shapes A, both for invariant coordinates and for generic 
coordinates. This generalizes known results for the rational normal curve A( n ), in particular 
the tight connection to eigenvectors of symmetric tensors; cf. [8j Ex. 5.2 and Cor. 8.7]. 

Section 6 discusses implications for the ED optimization problem (l2Tj) . The focus lies on 
locating all the critical points in Aa and in (Aa) v , and on certifying the global optimum. 
We work this out for several examples, including one from tensor decomposition over M. 


2 Duality for Binary Forms 


We fix a field K of characteristic zero and we represent univariate polynomials by binary 
forms. Let V = K[x,y\ n denote the space of binary forms of degree n and H v = K[u,v] n its 
dual vector space. For / £ V and g G H v , the pairing is defined by using partial derivatives: 


(f,9) 



8_ d_\ 
dx ’ dy ) 





d_ d_ 

du 1 dv 


g(u, v). 


(3) 


This is the scalar in K obtained by interpreting one polynomial as a differential operator 
and applying it to the other polynomial. Introducing coordinates on V and V v , we have: 


If / = E"=o G) a i x *y n 1 and 9 = E"=o 1 then (f,g) = E"=o ( 4 ) 

We regard / and g as elements in the projective spaces P(H) and P(H V ). Both of these 
spaces are identified with P n using homogeneous coordinates (a 0 : ■ ■ ■ : a n ) and (b 0 : • • • : b n ). 

A partition A of n is represented alternatively as a list of integers (Ai,..., A^) satisfying 
Ai > • • • > Xd > 1 and Y^t=i = n i or as a multiset {l mi ,... ,p mp } satisfying Xp=i j' m j = n - 
So, nij = #{i : A, = j}, and d is the number of parts of A, and p is the largest part of A. 
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The multiple root variety A a C P n comprises all binary forms / of degree n that have 
rrij roots of multiplicity j. Equivalently, writing 4 for linear forms, A a is the image of 


Bl ) d -g r\ (4,4,..., id) ^ / = 4 1 4 2 ---^ d - 


The variety A a has dimension d and degree as in ([ID. Among its smooth points are those 
given by d distinct linear forms 4- We determine the tangent space of A a at such a point. 

Lemma 2.1. The tangent space of A\ at a general smooth point f = Jlf=i equals 


t,a x = 


i =1 


h E P (K[x,y] d ) 


~ P d 


Proof. Write 4 = cqx + biy for some indeterminate (cq : bf) E P 1 . The tangent space TfA\ 
is spanned by the binary forms -£^f = xf jii and -J^f = y//4 for i = 1, 2,..., d. These 

generators have the required form if we take h(x,y) to be x-(J\* =l lj)/(-i or J/-(IT/=i4)/^- 

For the converse we note that K[x, y]d-i is spanned by {(IT/=i 4)/^» I i = 1, • • •, d). Indeed, 
this is a Lagrange basis, since the 4 are pairwise linearly independent. Hence, by multiplying 
these polynomials with x and y, we obtain a spanning set for the vector space K[x, y]d . □ 

Lemma 2.2. Let f be as above and g E P(L V ). Then g _L TfA\ if and only if 


d 



i =1 


d_ d_ 

du dv 


) 


annihilates g(u,v). 


(5) 


Proof. Let g(u,v ) be the binary form of degree d obtained from g(u,v) by applying the 
operator on the left of (j5p . Then g is zero if and only if g(u,v ) is annihilated by h(f^, 
for all h G K[x,y\ d if and only if g(u,v) is annihilated by {h ■ nf=i (Jjp J;) for all 
h G K[x,y\d- By Lemma [All this means that g(u,v ) is orthogonal to the space TfA\. □ 

The conormal variety of Aa, here denoted Coiia, is the Zariski closure of the set 
{(f,g) | / is a smooth point of Aa and g A TfA\} in P(V) x P(L V ) = P n x PT 


General results on projective duality ensure that ConA is an irreducible variety of dimension 
n— 1 in P n x P n . The dual variety (Aa) v is the image of the conormal variety ConA under 
projection onto the second factor P(L V ). This is an irreducible variety of dimension < n— 1. 
The Biduality Theorem um states that the conormal variety of (Aa) v coincides with Coiia, 
and hence ((Aa) v ) v = Aa- Lemma [2.21 implies the following description of the dual variety. 

Corollary 2.3. The points on the variety (Aa) v are binary forms g(u, v ) that are annihilated 
by some order n — d operator of the form nf=i where 4 ,...,^6 K[x, y]\. 

At this point, let us pause to illustrate the concepts seen above with a small example. 
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Example 2.4. Let n — 3 and take A = (3), the partition with a single part. The conormal 
variety Coii( 3 ) is the surface in P 2 x P 2 whose defining homogeneous prime ideal equals 

( a\ - aid 3 , aia 2 - a 0 a 3 , a\ - a 0 a 2 , 3 b\b 2 2 - 4 b 0 b% - 4 b\b z + 6 & 0 & 1 M 3 - /g', 

a o^o + 2+ o 2 & 2 , a 0 fei + 2ai& 2 + a 2 & 3 , cq&o T 2a 2 6i + o 3 & 2 , cq&i + 2a 2 6 2 + a 3 & 3 

The multiple root locus A( 3 ) is the twisted cubic curve in P 3 , consisting of cubes of linear 
forms, and defined by the first three quadrics in ([ 6 ]). Its dual is (A( 3 )) v = A( 2)1 ). It consists 
of binary cubics with a double root, and it is the discriminant surface in P 3 defined by the 
quartic in (JHJ). We have Coii( 3 ) = Con( 2 i i), after swapping the a-variables with the b- variables. 

We illustrate Corollary 12.31 for the point ((1 : 1 : 1 : 1), (1 : 2 : —5 : 8 )) in Con( 3 ). This 
represents the pair (f,g) where f — (x + y ) 3 and g = (u — v) 2 {u + Sv). Then g(u,v ) lies in 
(A( 3) ) v because (^+-^) 2 g(u, v ) = 0 , and f (x, y) lies in (A ( 2 j i)) v because (j%—§jj)f(x, y) = 0 . 
Note that the degree formulas ([I]) and ([ 2 ]) evaluate to 3 and 4 for A = (3). <0> 

The main result in this section is a parametric representation of the conormal variety. 

Theorem 2.5. The conormal variety Con A is the set of pairs ( f,g ) £ P n x P" of the form 

d d 

f(x,y) = Y[(Ux-Siy) Xi and g(u,v) = + Uv) n ~ Xi+2 ■ gi(u,v), (7) 

i= 1 <=l 

where (s* : tf) runs over P 1 and gi runs over binary forms of degree A* — 2. This parametriza- 
tion is a finite-to-one map (P 1 )" 11 x J\ — ■> Con A cP"x P" whose degree is mi!m 2 ! • • • m p \. 

In order for this theorem to make sense, we need to define the parameter space. We set 

J\ = Jou^P 1 x P Al_2 , P 1 x P a 2 ” 2 ) ..., P 1 x P Ad-2 ). 

This is the free join (or abstract join ) of the d — m\ varieties in the argument. Here 
is subtracted from d because the j-th argument is empty and gets deleted when A j = 1 . 
Note that J\ is the projective toric variety whose associated polytope is the free join of the 
product of simplices a i x <r\ 2 for j = 1,2 ,,d. The dimension of this join equals 

d 

dim(d7\) = (d — mi — 1) + max( 0,1 + (A* — 2 )) = n — m\ — 1 . 

2—1 

The point of the toric variety J\ is to ensure that (J7J) gives a well-defined rational map. We 
will see in Example 12. 101 that it is not a morphism. But that is not a problem since we can 
replace both the parameter space and image by their affine cones. The formulas in (J7J) are 
homogeneous. We use them in the next section to carry out the computations for Table [D 
By projection onto the second factor, Theorem 12.51 immediately yields a parametrization 
(P 1 )” 11 x J\ —-» (A a ) v C P n of the dual variety. Indeed, (A a ) v consists of all polynomials 
g(u,v) of the form in (IT]) . This can be rephrased in the language of projective geometry: 
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Corollary 2.6. The dual variety (Aa) v indexed by A is the join of d — mi multiple root loci 
determined by hook shapes, namely the loci n _ A . +2 \, where i — 1 ,..., d with A* > 2. 

Recall that a partition A is a hook if at most one part is not 1. This is special for us: 

Corollary 2.7. The dual variety (Aa) v is also a multiple root locus A M if and only if the 
partition A is a hook. In that case A = {l n ~ a , a} and p = {l a ~ 2 , n — a + 2} for some 
a G {2,... , n}. In particular, A a is self-dual whenever n is even and A = {l™/ 2 ^ 1 , n/ 2 + 1}. 

We easily derive the theorem from results that are well-known in commutative algebra. 

Proof of Theorem 12.51 We use apolarity as presented in the book of Iarrobino and Kanev 
[12]. Let / and g be binary forms, possibly of different degrees. We say that / is apolar to 
g if g[x,y ) is annihilated by the operator /(A ) A)_ Our Lemma [221 says that the tangent 

space TfA\ consists of all binary forms g of degree n such that is apolar 

to g. By m Lemma 1.31], this condition is equivalent to g having a generalized additive 
decomposition of the form g = ^ _ ( Al_1 ) +1 g l jQ~^ d ~^ +1 g d . Here, only terms with 

Aj > 2 may appear. This is precisely the representation on the right of ([7]), and we conclude 
that the proposed parametric representation of the conormal variety Coiia is correct. 

The parametrization is hnite-to-one because dim(CoiiA) — n — 1, which is the dimension 
of the parameter space (P 1 )" 11 x J\. Consider the fiber over a general point (f,g) in Con^. 
The first entry / is the image of precisely mi!m 2 ! • • -m p \ points in (P 1 )^. For each such point 
((si : ti),... ,(sd : td)), the second entry gives a homogeneous system of linear equations: 

d 

const. • g(u, v) = y '( s i u + tjv) n ~ Xi+2 ■ gj(u,v). (8) 

i=1 

The unknowns are the coefficients of gi,...,g d . The solution set is a linear subspace of 
J\. It is non-empty since (f,g) was assumed to lie in Con*. Since the parametrization is 
hnite-to-one, the solution space of the linear system must be a point. We conclude that the 
parametrization of the conormal variety Coiia given in (J7|) has degree • • -m p \. □ 

At this point it is important to note that our approach in this section is quite classical. 
All the ingredients are known and have been published elsewhere. What is new is their 
arrangement and interpretation. For instance, the parametrization (J7|) appears in [T2] but the 
connection to dual and conormal varieties was not made explicit. Likewise, our Lemma [2.11 
is among the statements in [T3], Theorem 7.1]. The dual variety does make an appearance in 
[Til Section 6] but it was not seen that the dual of a multiple loot locus for a hook is again 
such a locus. By putting all the known puzzle pieces together, we can go further in this 
paper. For instance, Katz states an inequality in [[Til Corollary 6.4], He conjectures in his 
introduction [Ti], p. 220] that equality holds. The following result proves Katz’ conjecture. 

Corollary 2.8. The variety of binary forms of degree n with a root of multiplicity a satisfies 

deg((A{in—a a }) v ) = deg(A{ 1 n-a+i a _ 1 }.) = (a — l)(n — a + 2). 
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Proof. The dual on the left is A{ 1 o-2 in _ a+2 } by Corollary 12.71 Using Hilbert’s formula ([I]), 
we see that both varieties have the same degree, namely (a — l)(n — a + 2). □ 

We close this section by recording the dimension of the dual variety, and by pointing out 
that our parametrization is in fact always identifiable , i.e. birational modulo permutations. 

Theorem 2.9. For any partition A of n, the dual variety (Aa) v has dimension n—mi — 1. 
The generalized additive decomposition in & represents a finite-to-one parametrization 
Jx —» (A a ) v C P n . This rational map has degree m 2 \ ■ ■ -m p \ and it is given by 

d 

^2(siU + tiv) n ~ Xi+2 ■ gfiu,v). (9) 

i= 1 
A^l 

Proof. It was shown in [IJJ Corollary 7.3] that dim((A,x) v ) = n — mi — 1. This is also the 
dimension of the parameter space J\. Theorem 12.51 implies that the map (jUJ) is finite-to-one. 

Fix a generic point g(u,v ) in the dual variety (A a ) v . Consider any smooth point f(x,y ) 
of the primal A^ at which g is tangent. The unordered set of linear forms ii that appear with 
multiplicity > 2 in f(x,y) can be recovered uniquely by Lemma [2721 Permuting linear forms 
that appear with the same multiplicity accounts for the size mfi. ■ ■ -m p \ of the fiber when 
writing out the multipliers (siU + tiv) n ~ Xi+2 in (jUJ). At this point, we still need to recover 
the forms (cq,..., gfi), but this can be done uniquely by the same argument as in the proof 
of Theorem 12.51 We conclude that the degree of the map ([9]) equals m 2 \ ■ ■ -m p \ □ 

Example 2.10. The rational parametrization in Theorem 12.91 can have base points, so it is 
generally not a morphism. Let A = (3,2), so (A^) v is the hypersurface in P 5 referred to as 
little apple in Example 14.21 The toric fourfold J\ is the free join of P 1 x P 1 and P 1 x P°. In 
these products, we Lx the points ((si : tfi), gfij and ((s 2 : t 2 ), gfij , where s± = s 2 = ti = t 2 = 1, 
gi = u + v, and g 2 = —1. We also Lx the scalar 1 for a point on the line that joins them. 
These choices specify a point in J\. Plugging into the formula §)]), we obtain 

g(u,v) = (siu + tiv) 4 ■ gi + (s 2 u + t 2 v) 5 ■ g 2 = (u + v) A -(u + v) + (u + v) 5 ■ (—1) = 0. 

Hence our point in J\ is a base point of the parametrization (M) for A = (3, 2). <0 

3 Equations, Multidegree, and More 

Given any parametrically represented variety, it is natural to ask for its implicitization. This 
concerns Lnding the ideals of polynomials that vanish on Con*, on A\, and on (Aa) v . For 
the multiple root loci, this is a well-studied problem [5] [15, 22]. Before reviewing what is 
known, we give the relevant dehnitions and we present a census of our varieties for n <7. 
We write J(Con^) for the ideal of the conormal variety in the N 2 -graded polynomial ring 

Q[a, b] = Q[a 0 , ai,..., a n , b 0 ,b i,..., b n \, with deg(eq) = (1, 0), deg(fei) = (0,1). 


^((f* : s *)Wi), ^ g(u,v) = 
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A 

Eqns of A a 

Multidegree 

Eqns of (A a ) v 

Hooks 

ED-degrees 

2 

2 

2,2 

2 

2 

2,4 

3 

2 3 

0,3,4 

4 

21 

3,7 

21 

4 

4,3,0 

2 3 

3 

3,7 

4 

2 6 

0 , 0 ,4, 6 

6 

211 

4,10 

31 

2 1 ,3 1 

0 , 6 , 6 , 0 

2 1 ,3 1 

31 

4,12 

211 

6 

6 ,4, 0,0 

2 6 

4 

4,10 

22 

3 7 

0,4, 6 , 3 

3 

4,4 

7,13 

5 

2 1U 

0,0, 0,5, 8 

8 

2111 

5,13 

41 

2 3 

0,0, 8 , 9,0 

46 

311 

5,17 

311 

4 6 

0,9, 8 , 0,0 

2 3 

41 

5,17 

2111 

8 

8 , 5, 0,0,0 

2 10 

5 

5,13 

221 

5 10 

0,12,16,6,0 

3 4 

5,5 

16,34 

32 

4 28 

0 , 0 , 12 , 21,12 

12 

41,5 

21,45 

6 

2 15 

0 , 0 , 0 , 0 , 6,10 

10 

21111 

6,16 

51 

2 6 

0 , 0 , 0 , 10 , 12,0 

4 1 ,5 3 ,6 1 

3111 

6,22 

411 

2 1 ,3 3 ,4 1 

0 , 0 , 12 , 12 , 0,0 

2 1 . S 3 ^ 1 

411 

6,24 

3111 

4 1 ,5 3 ,6 1 

0 , 12 , 10 , 0 , 0,0 

2 6 

51 

6,22 

21111 

10 

10 , 6 , 0 , 0 , 0,0 

2 15 

6 

6,16 

2211 

7 13 

0,24,30,10,0,0 

3 10 

6,6 

28,64 

222 

445 

0,0,8,16,12,4 

4 

6 , 6 , 6 

20,40 

33 

3 29 

0,0,0,9,18,12 

12 

51,51 

19,39 

321 

4 1 ,5 3 ,6 31 

0,0,36,56,24,0 

4 1 ,6 1 

51,6 

44,116 

42 

2 1 , 3 3 ,4 31 

0,0,0,16,30,18 

18 

411,6 

26,64 

7 

2 21 

0,0,0,0,0,7,12 

12 

211111 

7,19 

61 

2 10 

0,0,0,0,12,15,0 

6 10 

31111 

7,27 

511 

2 3 , 3 4 

0,0,0,15,16,0,0 

420 

4111 

7,31 

4111 

420 

0,0,16,15,0,0,0 

2 3 , 3 4 

511 

7,31 

31111 

6 10 

0,15,12,0,0,0,0 

2 10 

61 

7,27 

211111 

12 

12,7,0,0,0,0,0 

2 21 

7 

7,19 

22111 

9 16 

0,40,48,15,0,0,0 

3 20 

7,7 

43,103 

2221 

6 78 

0,0,32,60,40,10,0 

46 

7, 7,7 

62,142 

3211 

6 10 , 8 38 

0,0,72,105,40,0,0 

45, glO 

61,7 

73,217 

322 

g 364 

0,0,0,36,80,66,24 

24 

61,7,7 

94,206 

331 

3 10 

0,0,0,27,48,24,0 

7 8 

61,61 

39,99 

421 

4 20 , 6 42 

0,0,0,48,80,36,0 

8 1 , 10 3 , ll 2 , 12 49 

511,7 

52,164 

43 

3 10 , 4 66 

0,0,0,0,24,51,36 

36 

511,61 

51,111 

52 

2 3 , 3 4 ,4 38 

0,0,0,0,20,39,24 

24 

4111,7 

31,83 


Tabic 1: Multiple root loci, their duals, and their conormal varieties for n <7. 
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The ideal /(Con A ) is bihomogeneous and prime. According to textbook definition [T7[ §8.5], 
the multidegree of the N 2 -graded algebra Q[a, b]//(CoiiA) is a binary form of degree n + 1: 

Ca(^i, t 2 ) = + 62^2 1 + • • • + S n tft 2 - ( 10 ) 

The coefficients Si are nonnegative integers, known among geometers as the polar classes of 
the variety Aa. Geometrically, 5* is the number of intersection points in (L* x M{) n Coiia 
where L* is a generic plane of dimension i in P(I/) and Mi is a generic plane of codimension 
i— 1 in P(P V ). In particular, we have Si = 0 for i < codim(AA) and for i > dim((AA) v ) + 1. 

The ideal of the multiple root locus A a is obtained by eliminating the variables bo,... ,b n , 
and the ideal of its dual (A a ) v is obtained by eliminating the variables a 0 ,..., a n . In symbols, 

/(A a ) = /(Con A ) D Q[a] and /((A a ) v ) = /(Con A ) n Q[b], ( 11 ) 

These elimination ideals are N-graded and prime. Their degrees and codimensions in the 
respective P n can be read off from the first term and the last term of the multidegree: 

C x {t i,t 2 ) = deg(A A ) • C dim(AA) 4 im(AA)+1 + ••• + deg((A A ) v )-tf m((AA)V)+1 t c 2 ° dim((AA)V) . 

We computed the objects in (fTUl) and (ITIT) for all partitions up to n — 7. Our results 
are listed in Table [Q Each row corresponds to one partition A. The second column lists 
the degrees of the minimal generators of /(Aa). The third column lists the polar classes 
<5i, S 2 , ■ ■ ■, S n . The leftmost nonzero entry in that list is the degree of Aa and the rightmost 
nonzero entry is the degree of (Aa) v . The codimension of the variety in question is the 
number of consecutive left (resp. right) zeros plus one. The fourth column lists the degrees 
of the minimal generators of / ((Aa) v ) . The fifth column lists the corresponding collection of 
hooks {l Ai ~ 2 , n—Xi+ 2 } that make up the dual variety in the join construction of Corollarv l 2 . 6 l 

Finally, the last column refers to the Euclidean distance degrees in two coordinate sys¬ 
tems. This will be explained in Section 5. The second one is always Ca( 1, 1) = hi+h 2 + - ■ -+S n . 

We now discuss some general facts that a reader might discover by looking at Table [Q 
Whenever all hooks dual to A are identical then (Aa) v is a secant variety of that hook variety. 

Proposition 3.1. Suppose n > k{n — a+ 2). The k-th secant variety of the variety of binary 
forms of degree n with a root of multiplicity a is projectively dual to the multiple root locus 
with partition A = {i n - k ( n ~ a + 2 ) ^ (n — a + 2) k }. In symbols, 

&k (A{i»-«, a }) = (Aa) v . (12) 

This secant variety is non-defective: it has the expected dimension k{n — a + 1) + k — 1. 

Proof. This is the special case of Corollary 12.61 where A has the shape of a rectangle together 
with some singleton blocks, say A = {l n ~ ku ,u k }, and we set a = n — u + 2. The statement 
concerning the dimension follows from the first sentence in Theorem 12.91 □ 

The best known special case of this statement arises when a = n, or u = 2. The dual of 
Apn - 2 fc i 2 fc} is the /c-th secant variety of the rational normal curve A( n ). Its ideal is generated 
in degree k + 1, namely by the minors of a Hankel matrix. We see this in Table Q] for 
A = 21 , 211 , 22 , 2111 , 221 ,.... Also of interest is the case of a rectangular partition, when 
the secant variety (fT 2 ]i is a hypersurface. Oeding’s formula (J 2 ]) and Proposition 13.11 imply 








Corollary 3.2. The secant variety referred to in Propositio ns. 11 under the same hypothesis, 
is a hypersurface if and only if n = k(n — a + 2). The degree of that hypersurface is 

deg(a fc ( A {1 n-a >a} )) = (k + l)(n - a + l) fc . 

Example 3.3. Let k = 2, n = 6 , a = 5, so A = (3, 3). The variety (fliZD is a hypersurface of 
degree 12 in P 6 . It consists of all binary sextics / = £\£i -\-£\£± where the £i are linear forms. 
Consider also the case n — 12. Here the construction gives four interesting hypersurfaces: 


k 

a 

A 

forms 

degree 

2 

8 

( 6 , 6 ) 

^l 9 l + 

75 

3 

10 

(4,4,4) 

^l+^°<?2+4°<?3 

108 

4 

11 

(3, 3, 3, 3) 

Zfgi + ^2 92 + £fgz + ifg* 

80 

6 

12 

( 2 , 2 , 2 , 2 , 2 , 2 ) 

£? + £? + ■■■ + £f 

7 


Of course, the last hypersurface is defined by the determinant of a 7 x 7 Hankel matrix. <0 

The containment relation among multiple root loci is the order on partitions by refine¬ 
ment. Indeed, Aa C if and only if p refines A, i.e. every part of A is a sum of parts of p. 
Our next result characterizes the containment relation among their dual varieties. Note that 

(Aa) v C (A m ) v (13) 

cannot hold unless A has more l’s than p. for dimension reasons. Given a partition A = 
{l mi , 2 m2 ,... ,p mp }, we write A' for the partition {l m2 ,..., (p — l) mp }. Note that if A is a 
partition of n then A' is a partition of |A'| — n — d, where d = X mi is the number of parts. 

Proposition 3.4. Given two partitions A and p ofn, the inclusion fT3\) holds if and only if 
| A' | < |/i'| and, by adding to the parts, X' can be transformed to a partition X" refined by p!. 

Proof. This can be seen from Corollary 12.31 Let £ x ' be the differential operator that is 
used in that corollary to characterize (Aa') v - Here £ represents an arbitrary collection of 
linear forms. Our refinement condition means that every form g annihilated by £ x is also 
annihilated by £^ . This condition is equivalent to (JT3]) . More precisely, write m a for the 
number of parts in A', and for that in p'. If such a X" exists, then we claim that 

U Ann((4J A ') C U Ann((4j A ") C (J Ann((^ m J /i '). 

^m a £m a 

Here, is a /c-vector of linear forms and the exponent is written in multi-index notation. 
The first inclusion holds since each (£ ma ) x can be seen as a multiple of (£ ma ) x , and the 
second inclusion holds since being annihilated by (£ ma ) x can be seen as a special case of 
being annihilated by , where £ mi is obtained by duplicating some parts of £ ma . □ 

Example 3.5. We can check the inclusion (1X31) for various cases in Table [0 The inclusion 
holds for A = 2211 and p = 321. Indeed, A' = 11, and p' = 21 refines X" = 21 . Likewise, 
if A = 321 and p = 222 then p' = 111 refines X! = A" = 21. This explains the inclusions 
(A 22 ii) v C (A 32 i) v C (A 222 ) v C P 6 . It follows that the unique quartic polynomial 4 1 that 
vanishes on (A 3 2 i) v must be the 4x4 Hankel determinant whose hypersurface is (A 222 ) v - § 
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We now discuss some of the entries in the column Eqns of A a in Table [T] and thereby 
review the literature on the ideals /(A a ). For A — (n), we see the ("’) quadrics that define 
the rational number curve, and for A = {l n ~ 2 , 2} we get the discriminant of degree 2n —2. If 
a > |n/2j + 2 then /(A^n-a^}) is generated in degree < 4, as shown by Weyrnan [22] . The 
case a = |_ri/ 2 j +1 is of special interest: here is the nullcone, i.e. the variety defined 

set-theoretically by all SL 2 -invariants of binary n-ics. Note that the nullcone is self-dual 
when n is even (by Corollary 12 .7|h and its ideal is generated by quartics when n is odd. 

Chipalkatti [5J Conjecture 6.1] had conjectured that /(A a ) is generated in degree < 4 
whenever A has only d = 2 parts, and this was proved by Abdesselam and Chipalkatti in 
|2J Proposition 20]. For further details on their approach see jl] Section 7]. The papers 
d 121 El |22] stress the fact that all our ideals are invariant under the action of SL 2 , and hence 
each of their graded components are direct sums of irreducible SL 2 -modules. Chipalkatti 
describes the minimal free resolutions of J(Aa) in terms of SL 2 -modules for A = (3,2) in [3 
§4.1] and for A = (3,3) in [0 §4.2], He discusses the 364 sextics for A = (3,3,2) in [5] §3.1]. 

The ideal /((Aa) v ) of the dual variety has been studied only in two cases, namely when 
A is a hook (by self-duality) and when A has only parts 1 and 2. In the latter case, it is 
generated by minors of a Hankel matrix. In all other cases, little seems to be known. Of 
course, /((Aa) v ) will often be principal (when m i = 0), and we seek to find the generator 
of such an ideal explicitly. This is accomplished for some interesting cases in Section 4. 


4 Real Rank Boundaries 

In this section we present an application of our duality theory to the study of real ranks 
of binary forms PIE]. Let / G M[x, y\ n be a general binary form of degree n with real 
coefficients. The complex rank of / equals r = |~(n + l)/2]. This means that / is a sum of r 
powers of linear forms over C, but not fewer. We consider the set 7 Z n of all real binary forms 
/ that admit a rank r decomposition also over R. In other words, 7 Z n is the set of all forms 
/ G R[x, y\ n such that both the real rank of / and the complex rank of / are equal to r. 

The set 7 Z n is a full-dimensional semi-algebraic set inside R[x, y] n . Its topological boundary 
dIZ n is the set-theoretic difference of the closure of 7 Z n minus the interior of the closure of 
7 Z n . Thus, if / G d7Z n then every open neighborhood of / contains a generic form of real 
rank equal to r and also a generic form of real rank bigger than r. The topological boundary 
<97 Z n is a semi-algebraic subset of pure codimension one inside the real affine space R[x, y] n . 

We define the real rank boundary , denoted <9 alg (7 Z n ), to be the Zariski closure of the 
topological boundary <97 Z n . We view c7 a ig(7^ n ) as a closed subvariety in the complex projective 
space P(C[x, y ] n ) = R n - It is pure codimension one, so it is defined by a unique (up to scaling) 
squarefree polynomial in the coordinates ao, ai,..., a n on P n . We compute that polynomial: 

Theorem 4.1. Let n > 5. If n = 2k — 1 is odd then the real rank boundary d a \ g (fJZ n ) is an 
irreducible hypersurface of degree 2 k{k — 1) in P n . This hypersurface is the dual (Aa) v of the 
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multiple root locus A a where A = {2 k 2 , 3}. Its defining polynomial is the discriminant of 


q(u,v ) 


det 


u k 

u L v 

1 

• • • uv 

v k 

Oq 

ai 

Ofc-l 


a\ 

a 2 


Ofc+1 

a 2 

a 3 

Ofc+1 

Ofc+2 

Ofc-l 


1 



(14) 


If n = 2k is even then the real rank boundary d a \ g (IZ n ) is the union of the two irreducible 
hypersurfaces (Aa) v where A is { 2 fc ~ 3 , 3 2 } or { 2 fc ~ 2 ,4}. Their degrees are 2k(k — 1 )(k — 2 ) 
and 3 k[k — 1). These arise as irreducible factors of the Hurwitz form of the discriminant 
Ajifc- 1 , 2 } when evaluated at the line in P(M[x, y]k+i) that is the kernel of the kx(k+2)-matrix 



cq 

eq 



a 3 

a 2 

‘ < 1 k+ 1 

&k+ 2 

(15) 



®n— 1 

) 



A third irreducible factor appearing in this specialized Hurwitz form is the Hankel determinant 
(A{ 2 fcj) v , which has degree k + 1, but this is not part of the real rank boundary d a \ g {IZ n ). 

For the relevant background on Chow forms and Hurwitz forms we refer to m and [21] • 
Before presenting the proof of Theorem 14.11 here is an illustration of the first three cases. 


Example 4.2 (Little Apple and Big Apple). Let n = 5, so k = 3 and A = (3, 2). In this first 
case, Theorem 14.11 was proved by Comon and Ottaviani in [FJ §5]. The polynomial defining 
the hypersurface (Aa) v C P 5 is their apple invariant Ji 2 . It has 228 terms of degree 12. 

The next odd case is n = 7, so k = 4 and A = (3,2,2). The real rank boundary (Aa) v 
is a hypersurface of degree 24 in P', namely the join of the surface A( 6>1 ) and the threefold 
<t 2 (A( 7 )) = er 2 (i/ 7 (P 1 )). Its defining polynomial is computed via (fT4]b It has 38082 terms. <0 

Example 4.3 (Chow and Hurwitz forms). Let n = 6 . We consider the discriminant A( 211 ) 
of binary forms of degree k + 1 = 4. This discriminant is a threefold of degree 6 in P 1 . 
Its singular locus consists of the surfaces A( 3;1 ) and A( 2)2 p The Hurwitz form of A( 2j11 ) is a 
polynomial of degree 30 = 3 -6 + 2-4 + 4 in the ten Pliicker coordinates on the Grassmannian 
of lines in P 1 . It is reducible because A( 2 i i i i) is singular in codimension 1. It factors as 

Hur(A( 2 ; i i i)) = Chow(A( 3 i i )) 3 • Chow(A ( 2]2 )) 2 • Tan. (16) 

Here Cliow(X) denotes the Chow form of a surface X in P 4 , i.e. the polynomial in Pliicker 
coordinates that vanishes when the line intersects X. The degree of the Chow form equals 
the degree of X, so it is 6 and 4 in our two cases. The last factor Tan is the condition for the 
line to be tangent at a smooth point of A( 2j1;1 ). This is a quartic in the Pliicker coordinates. 
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We now take our line in P 4 to be the kernel of the 3x5 matrix in (TT5T) . More precisely, in 
the formula (fTUj) we substitute the 10 maximal minors of (1131) for the 10 Pliicker coordinates, 
ffere signs have to be taken into consideration carefully. Each maximal minor is a cubic, so 
the degrees above have to be tripled. Equation (TT 6 |) now has degree 90 = 3-18 + 2-12 + 3- 4. 

The polynomial obtained from Chow(A( 3]1 )) is (A( 4]2 )) v . It has 3140 terms of degree 18. 
The polynomial obtained from Chow(A( 2 ; 2 )) equals (A( 3]3 )) v . It has 560 terms of degree 12 . 
The degree 12 polynomial obtained from Tan is the third power of the Hankel determinant 


(A( 2 , 2 , 2 )) v — det 


fa 0 

a.\ 

02 

03 \ 

a\ 

02 

03 

a 4 

02 

03 

a 4 

05 

\03 

a 4 

a,5 

aj 


(17) 


To illustrate the statement of Theorem 14.11 we present two explicit forms that lie in dTZ$. 

The binary form f = y 6 + 15x 4 y 2 lies in <977 6 D (A( 4j2 )) v but not in (A( 3j3 )) v U (A( 2 i 2 , 2 )) v - 
Its line of apolar quartics is Lf = {s-wu 3 + t ■ {u 4 — v 4 ) | (s : t) G P 1 }. This has discriminant 
(27s 4 + 256t 4 )t 2 , with only real root at (s : t) — (1 : 0). That quartic lies in k = A( 3j1 ), and 
it is a limit of quartics with four real roots, and also a limit of quartics with two real roots. 
From this we can construct generic sextics f± e close to / whose real ranks are 4 and 5. 

The binary form / = y 6 + 5x 2 y 4 — 5x 4 y 2 — x 6 lies in <977 6 n(A( 3!3 )) v but not in (A (4i2 )) v U 
(A( 2 , 2 , 2 )) v - Its line of apolar quartics is Lf = {s-(w— v) 2 (u+v) 2 + t-uv{u 2 +v 2 ) | (s : t) G P 1 }. 
This has discriminant (16s 2 + t 2 ) 2 t 2 , with only real root (s : t) — (1 : 0). That quartic lies 
in v = A( 2 , 2 ) • We can construct nearby generic sextics f± e whose real ranks are 4 and 5. 

The proof below will explain the relevance of Lf, k, and v for <977-6. It will also show why 
the Hankel determinant (A( 2 , 2 , 2 )) V meets the boundary of 1Z§ only in lower dimension. <0 

Remark 4.4. It is instructive to see the factorization (1T61) using iterative discriminants. Let 
A(x) and B(x) be univariate quartics, and consider A(x)+tB(x) where t is an unknown. Then 


discr t (discr x (A(x) + tB(x))) = ( [6, 6] AjB ) 3 • ([4,4] AB ) 2 - ([4,4 ]a,b), 


where [e, e] a,b stands for a big expression that has degree e in the coefficients of A and of B. 


Proof of Theorem f.l 
case n = 2k 


We write / for a generic form in M[x,|/] r 


We first consider the odd 

1. The apolar ideal, consisting of all forms in M[w,n] such that 
annihilates f(x,y), is generated by forms q(u,v ) of degree k and r(u,v ) of degree k + 1. 
Furthermore, the dual form q(u,v) has the Hankel determinantal representation in (fT4]b 
Since / is generic, there exists a unique decomposition f(x,y) = J^^ = 1 (si^ + tiy) n . The 
points (sj : tf) are the complex roots of q. Hence / lies in 7 Z n precisely when q is real-rooted. 
By real-rooted we mean that q is square-free and its roots are all real. 

Suppose now that the form / moves and passes through the boundary of 7 Z n . Then two 
real roots of q(u, v ) merge and become a double root. At this point, the discriminant of q(u, v ) 
vanishes. Corollary 12.31 implies that this discriminant equals (A^) v where A = {2 fc_ 2 ,3}. 

Next consider the even case n = 2k. The apolar ideal of / is generated by two forms of 
degree k + 1. Let Lf be the line in P fc+1 spanned by these two forms. The points q on Lf 
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correspond to the distinct decompositions f(x,y) = Y2i=i( s i x + Uy) n - Namely, the (sj : t,) 
are the roots of q. This means that / lies in 7 Z n if and only if some q in Lf is real-rooted. 

The set of all real-rooted forms is a connected full-dimensional semi-algebraic subset 
of P^ +1 , and its algebraic boundary is the discriminant 5 = Apun. Note that 5 is a 
hypersurface of degree 2k. Its singular locus consists of the codimension two loci k = A^fc -2 3 } 
and v = 3 2 2 }. Their degrees are deg(«:) = 3 (k — 1) and deg(V) = 2[k — 1 )(k — 2). 

Suppose that the form / moves along a general curve. Consider its image under the 
map / 1 —Lf into the Grassmannian of lines in P^ +1 , here denoted Gr. As / crosses the 
boundary of 7 Z n , the line Lf transitions from intersecting to not intersecting the subset of 
real-rooted forms. At the transition point, the form / is in <9 a i g (7 Z n ). One of the following 
three scenarios might happen: (i) Lf intersects k, or (ii) Lf intersects is, or (iii) Lf is tangent 
to 5 at a smooth point. Each of these conditions defines an irreducible hypersurface in the 
Grassmannian Gr. Their equations are the irreducible factors in the Hurwitz form Hur(<5). 
The Hurwitz form of the discriminant factors as follows in the coordinate ring of Gr: 

Hur(5) = ChowtA ) 3 • Chow(V ) 2 • Tan(5). (18) 

The exponents 3 and 2 arise from the classical Pliicker formula for the dual of a plane curve. 
Now, since 6 is a hypersurface of degree 2k, its Hurwitz form has degree 2k(2k — 1). Hence 

2k{2k— 1) = 3deg(ft) + 2&eg(is) + deg(Tan(<5)) = 9(/c — 1) + 4(/c — 1) (/c — 2) + deg(Tan(<5)). 

This uses the fact that the degree of a Chow form in Pliicker coordinates equals the degree of 
its variety. We conclude that Tan(h) is a polynomial of degree k + 1 in Pliicker coordinates. 

The map P n ---» Gr, / h -> Lf is birational. Indeed, every generic line has the form Lf for 
a unique sextic / that is recovered by solving the differential equations represented by Lf. 
Moreover, the base locus of this inverse map is precisely the tangential hypersurface Tan(5). 

We now examine what happens to the three irreducible factors in (fT8l) when pulled back 
under the map P" —Gr. Algebraically, the line Lf is given as the kernel of the matrix 
m So, to compute the pullback of (1181b we need to replace the Pliicker coordinates by the 
(appropriately signed and scaled) maximal minors of (fT51h These minors are polynomials of 
degree k, so we must multiply the above degrees by k. This gives the degrees 3 k{k — 1), 
2 k(k — 1 )(k — 2 ) and k{k + 1 ) for the pullbacks of the three irreducible factors in (fT 8 ]h 
The first two pullbacks are irreducible as polynomials in R[a] = M[ao,..., a n \. They are 

Chow(Ac)(L/) = (A{ 2 fc- 2 j4 }) v and Chow(zz)(L/) = (A{ 2 fc-3 i3 2j) v . (19) 

As before, this follows from the description of the dual hypersurfaces in Corollary 12.31 
The third factor in (TT8|) becomes a reducible polynomial in M[a], Namely, we have 

Tati (L,) = ((A, 2n ) v )t (20) 

In words: the pullback of Tan is the kt\i power of the Hankel determinant of order k + 1. 

This can be seen as follows. A general point / on the Hankel hypersurface satisfies 
rankc(/) = k. It is annihilated by some binary form g(u,v) of degree k, and the line Lf 
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is spanned by u ■ g(u,v ) and v ■ g(u,v). This means that Lf is tangent to S at k points 
over C, given by the k roots of g(u,v). We see that the pullback Tan (Lf) contains (A{ 2 fe}) v 
set-theoretically. Both are irreducible varieties, and hence they are equal as sets in P n . By 
comparing degrees, we conclude that f[20j) holds. 

Our argument also shows the following fact: if a line of the form Lf is tangent to the 
discriminant 6 at one smooth point then it is tangent to <5 at a scheme of length k. 

We have proved that <9 a i g (7£ n ) has at most three irreducible factors. However, the correct 
number is two. The two hypersurfaces that appear in <9 a i g (7 Z n ) are those in dTU]l . This was 
shown already for n = 6 . Towards the end of Example 14.31 we exhibited one binary form / 
for each of these two boundary strata. In general, a construction can be made as follows. 

We consider lines in the space P fc+1 of binary forms of degree k + 1 that look like 

L k = {s • uv 3 f(u,v) + t ■ (u 2 + v 2 )g(u,v) \ (s : t) G P 1 } and 

L v = {s • (u — v) 2 (u + v) 2 f{u , v) + t ■ (u 2 + v 2 )g{u , v) \ (s : t) G P 1 }, 

where / and g are generic real-rooted binary forms of degrees k — 3 and k — 1 respectively. 
The line L K meets the discriminant S only in the cusp locus, and the line L v meets S only 
in the node locus. The binary form has k — 1 distinct real roots at these intersection points. 
Now, for suitable choices of / and g, the discriminant of the pencil of binary forms has no 
real roots other than (s : t) = (1 : 0). We fix such / and g, and we set L = L K or L = L v . 

The pencil L consists of binary forms of degree k + 1 that have precisely k — 1 real roots. 
Moreover, L is not contained in the hypersurface Tan. Recall that the inverse to P" --■> Gr 
is well-defined in a neighborhood U of L. Hence, for each U EU there exists a unique form 
fu of degree n such that Lf v = U. By construction, there exist generic points U± e in U such 
that U +e contains a real-rooted form, and U- e contains no real-rooted form. Then f +e = fu +e 
has real rank k + 1 while /_ e = fu_ t has real rank > k + 2. We conclude that / is in d1Z n . 

It remains to be seen that the Hankel determinant is not a factor of the real rank 

boundary <9 a i g (7 Z n ). The proof is by contradiction. Suppose that / G <97 Z n has rank c (/) = k, 
and write / = JA =1 The ii are linear forms over C, and these are unique. We can 
approximate / by a sequence of forms in 7 Z n . Each is a sum of k + 1 powers of real linear 
forms. A convergence argument shows that A, • • •, Ac must have real coefficients as well. 

Fix a generic real linear form f 0 and consider f + e £q and / — where e > 0 is small. The 
square Hankel matrices corresponding to these two forms are invertible. Their determinants 
have opposite signs. This can be seen from the Matrix Determinant Lemma. The two forms 
lie on different sides of the Hankel hypersurface. Both have real rank k + 1, and hence both 
lie in 7Z n . This contradicts our assumption that / is in the topological boundary of 7 Z n . □ 


Remark 4.5. The real projective space P^ of binary forms of degree n is stratified by real 
rank. Let r = \[n +1)/2]. Blekherman [3] showed that each integer in {r, r +1,..., n — 1, n} 
arises as the real rank of some open stratum. It would be very interesting to determine the 
algebraic boundary that separates real rank i from real rank i + 1 , for any i G {r,..., n — 1 }. 
Currently, only the two extreme cases are known. The algebraic boundary for i — n — 1 is 
the discriminant A { ln -2 2 }, by [ 6 j Prop. 3.1], and the case i = r is resolved by Theorem 14.11 


14 





5 Euclidean Distance Degrees 

This section is concerned with the following optimization problem over the real field K = R. 
Let h G R[x,y] n be a fixed binary form. We seek / e A a that is closest to h. In symbols: 

minimize (/ — h, f — h) subject to / e A\. (21) 

We call this the special ED problem when the inner product is as in (j3J) and (j4j). Alternatively, 
we may also use a positive definite quadratic form that is generic. The resulting scenario 
is the generic ED problem for A\. For instance, if we replace (j3J) with (f,g) = Y^=o 7i°A, 
where 70 ,7i, ■ ■ ■, 7n are random positive reals, then this leads to a generic ED problem. 

The map that takes the given input h to the optimal solution f* = f*(h ) of the problem 
m is an algebraic function. The degree of that algebraic function is known as the ED degree 
of the variety A*. For an introduction to this topic and its basic results see PJ. A follow-up 
study, aimed at varieties that admit a determinantal representation, was undertaken in [20] . 

In the algebraic approach to solving (12T| one writes the critical equations using Lagrange 
multipliers and one removes the singular locus. We shall see concrete examples for Aa 
shortly. The ED degree is the number of complex solutions to these equations for generic 
data h. Of course, the optimal point f* is real and it is among these complex critical points. 

Varieties such as Aa come with a natural invariant coordinate system, and our special 
inner product (l3l)-(]TJ) gives the standard Euclidean distance for that coordinate system. The 
corresponding ED degree is the special ED degree of Aa- By contrast, the generic ED degree 
of Aa is usually larger. This is the degree for the generic ED problem described above. 

In the last column of Table [Q we list the special ED degree and the generic ED degree. 
The first pair in each box, corresponding to the rational normal curve A( n ), equals n, 3n — 2. 
These numbers were derived, for arbitrary n, in [SJ Example 5.12] and P Corollary 8.7]. 
Both notions of ED degrees are preserved under duality [SJ Theorem 5.2], and the generic 
ED degree coincides with the sum of the polar classes, C\(l, 1) = E A by pi, Theorem 5.4], 
Some of the ED degrees in Table |T] appear in [20]. For instance, the ED degrees 7,13 
for A = (2,2) appear in [2UJ Example 1.1]. When A is dual to a collections of hooks of 
the form n,n or n, n,n, the variety (Aa) v is given by Hankel matrices of rank 2 or 3. The 
corresponding ED degrees in Table |T] are found in [2Dj Table 4] . Therein, the left table for 
A = Q n gives generic ED degree, while the right table for A = O n gives special ED degree. 

From Table [T] we can guess the two ED degrees for the variety of binary forms of degree 
n that have a root of multiplicity a. Our next goal is to prove that this guess is correct. 

Theorem 5.1. For any hook A = {l n ~ a , a], the special ED degree of the variety Aa is always 
n, independently of a, whereas the generic ED degree of A a equals (2 a — 1 )n — 2(a — l) 2 . 

Proof. We begin with the generic ED degree. It is the sum of the polar classes 8^. The 
codimension a — 1 of Ajp-o^j equals the dimension of its dual (A^n-a a }) v = A^ia- 2 n _ a+2 y. 
This means that precisely two polar classes are nonzero, and the generic ED degree equals 

8 a—i T 8 a = deg(A{ 1 n-a ia j) + deg(A{ 1 a -2 n _ a+2 }) = a(n — a + 1) + (n — a + 2)(a — 1), (22) 

by Hilbert’s formula ([Tj) . This expression equals (2a — 1 )n — 2(a — l) 2 , as desired. 


15 


Next consider the special ED degree. Our goal is to compute the complex critical points 
f* of the optimization problem (1211) where h is a fixed generic binary form in V = M[x, y\ n . 
We now identify V with its dual space V v by way of the distinguished inner product ((2D, and 
we regard the conormal variety Con A as an affine cone of dimension n+ 1 in V x V = R 2n+2 . 
By ED duality [HI §5], our problem is equivalent to solving the system of linear equations 

/ + g = h for (/, g) G Con A . (23) 

These n + 1 inhomogeneous linear equations have finitely many complex solutions ( f*,g*) 
on the (n + l)-dimensional affine variety Con A . Their number is the special ED degree. 

We now write the equations (1231) using the parametrization of Con A given in Theorem 12.51 

h(x,y) = (tx - sy) a ■ gi(x,y) + (sx + ty) n ~ a+2 ■ g 2 (x, y). (24) 


The two summand are unknown points in Ann-a^} and in A{!a- 2 iri _ a+ 2 j. Both are now 
regarded as affine varieties in V = C n+1 . The forms g\ and g 2 are unknown and they have 
degrees n — a and a — 2 respectively. Furthermore, (s : t) is an unknown point in P 1 . Hence 
there are n + 1 unknown parameters in total, to match the n + 1 given coefficients of h(x, y ). 
Thus (124D is a square polynomial system in C n+1 , and we need to count its solutions. 

Any representation (1241) over C of the given binary form h(x, y) will be called an orthogo¬ 
nal decomposition of type a. Thus, our proof reduces to establishing the following assertion: 
a general binary form of degree n has precisely n orthogonal decompositions of type a. 

In what follows we describe the binary form whose zeros are the points (s : t) that can 
occur in (|24|) . In other words, we eliminate the unknown binary forms g\ and g-z from (124)) . 
Once (s : t ) is known, the coefficients of g\ and g 2 can be recovered by solving a linear system 
of equations. So, it suffices to identify that binary form and to show it has degree n. 

We define an endomorphism L^ of the vector space V = R[oi, y] n as follows: 


h ik) m*,y) = 

i =0 


U\ f) k 


This linear differential operator generalizes the S 02 -invariant vector held 


= 
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d_ 

dy 
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d_ 

dx 


and the second order operator 
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n{n — 1 ) [ dy 
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2 xy 


d 2 

dxdy 


+ y' 


dx 2 


The linear map L^ : V —> V can be written explicitly in terms of coordinates as follows. 
If f = Y^i=o {T) aiXl y n ~ l then the coefficients of Lf k \f) = Xq =0 (^)bjX j y n ~ j are 



( 25 ) 
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( 26 ) 


If we apply the n-th order operator then this amounts to a rotation by 90 degrees: 

(L<”>(/))(z,y) = f(-y,x). 

Now, the relevance of the differential operator L^ for our proof is as follows: 

Lemma 5.2. Let h be a binary form of degree n. Suppose that (s : t) G P 1 occurs in an 
orthogonal decomposition p4 \ ) of type a. Then (s : t) is a root of the binary form ZA° -1 )(/i). 

Lemma 15.21 will be proved further below. We first derive the theorem from the lemma. 
The linear map ZA a_1 ) is not zero. Hence, for generic h, the binary form l/ a_1 )(h) is 
nonzero and has degree n. Such a binary form has at most n distinct roots. Therefore, 
by Lemma 15.21 we know that the special ED degree is at most n. What we must prove is 
that the special ED degree is at least n. We do this by exhibiting, for every n and a, one 
particular binary form of degree n that has n distinct orthogonal decompositions of type a. 
We begin with the observation that the following two binary forms have real coefficients: 

K(x,y) = \{(x + V^i-y) n + {x- V^i-y) n ), (27) 

k n ( x ,y) = 9 J-^ ((a + -y) n + {x- \/ z I-2/) n ). 

These two forms are invariant under rotation by 2n/n, i.e. they are fixed by the endomor¬ 
phism p : V —> V that maps x to (cos —)x — (sin —)y and y to (sin —)x + (cos —)y. 
Case 1. Suppose that n is odd. For our special binary form we take 

hn(x,y) = x 11 - (” )x n ~ 2 y 2 + (l)x n ~ A y A -± (™_^xy n ~ l . 

Note that all the exponents of x are odd while that of y are even. For odd a we have 

h = x a gi + 9 ”-“ +2 92 . (28) 

for suitable forms g\ and g 2 depending on a. Similarly, for even a we have 

h = y a gi + x n ~ a+2 g 2 . 

For each decomposition, we obtain n — 1 others by acting with the rotations p, p 2 ,, p n ~ x . 
Case 2. Suppose that n even and a is even. We also take 

h n (x, y) = x n — (™)x n ~ 2 y 2 + { n ^x n ~ A y A - ±y n . 

Then (|28l) also holds. Now, h n is fixed by p n//2 , so applying the rotations p 1 gives only n/2 
distinct orthogonal decompositions of type a. However, we have h n (x,y) = ±h n (y,x), so by 
permuting the two variables we get additional decompositions. These are not equal to any 
of the previous ones. In total, this yields ^ • 2 — n distinct decompositions (l24f) for h n . 
Case 3. Suppose that n even and a is odd. In that case we take 

k n (x, y) = i^)x n - 1 y - Q)^-y + Qy-y - ■ ■ ■ ± { n \)xy n -\ 

and the argument is the same as in Case 2. This completes the proof of Theorem 15.11 □ 
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We now turn to Lemma [5721 Note that this result is familiar in the case a = 2 , when (l 2 Tj) 
asks for the best rank 1 approximation of a symmetric 2x2x • ■ ■ x2 tensor h. Finding that 
approximation amounts to computing the eigenvectors of h\ see [HI Corollary 8.7]. However, 

by definition, the eigenvectors of h are the roots of L^\h) = (1 /n) ■ det (q^q x QhjQy 

Proof of Lemma\5M We claim that the operator ZW satisfies the following identity 

L (k) {f(x,y)) = (-l)- k . L ^- k \f(-y,x)) for k = 0,1,..., n. (29) 


The special case k = n is (1261) . One can show that (129|) holds by a direct computation, using 
the formula (1251) . Equivalently, we check that it holds for monomials and extend by linearity. 

Suppose now that (s : t) G P 1 occurs in an orthogonal decomposition (j24l) of type a. We 
apply the differential operator L (“ _1 ) to both sides of that equation. This implies 

L {a ' l) (h(x,y)) = L (a_1) {{tx - sy) a gi(x,y)) + L (a_1) ( (sx + ty) n ~ a+2 g 2 (x, y )). 

Applying (1291) to the summand on the right, we conclude 

U"- 1 ) (/>(*,»)) = V‘- 1 \(tx-syr gi (x,y)) + (—1 )“-“+ 1 L<"-“+ 1 )((ii-s !/ )"-“+ 2 92 (-9.i)). 

In both summands, a k -th order differential operator L^ is applied to a binary form that 
has (s : t) as a root of multiplicity at least k + 1. Both of the resulting binary forms still 
have (s : t ) among their roots. This shows that iA a_ b (h(x, y)) vanishes at (s : t). □ 


The difference between the two ED degrees in Theorem 15.11 is 2(n — a + l)(a — 1). We 
shall explain this number and how to think about the operator lA a_1 h If we fix values for the 
parameters s and f, then the equation ( 1 ^I|) translates into an inhomogeneous linear system 
of equations whose unknowns are the coefficients of gi and g 2 • We have n+ 1 linear equations 
in n — (n — a + 1) + (a — 1) unknowns. There is no solution for generic s and t. 

We write our inhomogeneous linear system as an {n + 1) x (n + l)-matrix M. n , a - The 
first row consists of the coefficients of h. The next n — a + 1 rows contain the monomials of 
degree a in (s, t), and the last a — 1 rows contain the monomials of degree n — a + 2 in (s, t). 
We seek row vectors of the form (—l,^, g 2 ) that he in the left kernel of M. n , a - The matrix 
has a banded structure, like Sylvester’s matrix for the resultant. For instance, for A = (3,1), 


AZ 4,3 — 


In order for our linear equations to have a solution, the determinant of M. n>a must be zero. 

The degree of det(Af nia ) in (s,t) is precisely the generic ED degree. This is best seen 
from (122|) . If we replace ([3]) by a generic inner product then the rows for g\ change by a linear 


( ho 

h 

h 2 

h 3 

K\ 

—s 3 

3 s 2 t 

—3 st 2 

t 3 

0 

0 

-s 3 

3 s 2 t 

—3 st 2 

t 3 

t 3 

3 st 1 

3 s 2 t 

s 3 

0 

^ 0 

t 3 

3 st 2 

3 s 2 t 

s 3 J 


18 








transformation, whilst the rows for g 2 change by the inverse linear transformation. Thus, for 
the generic ED problem, the critical points are precisely the roots of the determinant. 

ffowever, for the special ED problem, our determinant admits the following factorization: 

det (M n: a) = ±(L^ a - 1 \h))(s,t)-(s 2 + t 2 )^- a+1 ^ a - 1 l (30) 

The binary form s 2 + t 2 vanishes if and only if the last n rows of M. n ,a are linearly dependent. 
The degree of the extraneous factor is 2(n — a + l)(a — 1), and subtracting that from the 
generic ED degree gives n. The remaining factor is the remarkable binary form Z,(° - 1 )(/j). 

We close this section with a conjecture, namely that the two ED degrees of A a always 
have the same parity. According to Table [H this holds for all partitions with n < 7. For 
hooks, it is proved by Theorem 15.11 and the underlying reason is seen clearly in (1301) . In 
general, the extraneous components should come from isotropic quadrics like s 2 + f 2 , and 
parallelities like sv — tu. These quadrics suggest that the difference in ED degrees is even for 
all A. This conjecture is related to [201 eqn. (3.5)]. At present we do not know how prove it. 


6 ED Duality in Action 

We now illustrate how our results can be applied to End exact solutions to the optimization 
problem (121]) . Following [ 8 ] [20], our approach is to compute all critical points and then select 
the best real critical point. By ED duality [ 8 J §5], the critical points are found by solving 
linear equations on the conormal variety. Given h, we need to compute all decompositions 

h(x,y) = f(x,y) + g{x,y) where (f,g)e Con A . (31) 

If h is generic then the number of such decompositions is the special ED degree. The binary 
forms / that arise in the decompositions (I3T1) are precisely the critical points on Aa for the 
Euclidean distance to h, and similarly the forms g are the critical points on its dual (Aa) v . 
The proximity of the solution is reversed under duality because \\h — f\\ 2 + \ \h — g\\ 2 = \\h\\ 2 . 
This follows from h = f + g and (f,g) = 0. In particular, if / is the closest point to h among 
those on Aa, then it is paired with the farthest critical point g on (Aa) v , and vice versa. 

In what follows we illustrate how one might solve (l3Tj) and hence (I2T]) in practice. We 
begin with n — 5 and A = (3,1,1). For our given data point we take the binary quintic 

h(x,y) = ^{x + V-l-y) 5 + 7 j{x~ \/-lw) 5 + y 5 = x 5 - 10x 3 y 2 + hxy 4 + y 5 . (32) 

This is a slight variant of ([27]) . The primal problem is to hnd the closest quintic / G A( 3 ) i ! i) 
with a triple root, and the dual problem is to find the closest quintic g G A( 4ii ) with a 
quadruple root. By Theorem 15.11 the equation ([3lT) has five solutions on Coiia- They are 
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Nk 

f = h-g 

9 = h- f 

i mk 

3.02 

(1.471a; 2 - 5.582xy + 3.585 y 2 )(x + 0.785 y)' 6 

(-1.238x - 0.735y)(0.785x - yf 

13.98 

4.92 

(-0.263x 2 - 0.211 xy - 0.338 y 2 ){x - 3.132y) 3 

(0.0131x - 0.403y)(3.132a? + y) 4 

12.08 

5.09 

(-0.263x 2 + 0.167 xy + 0.0346y 2 )(x + 3.020y) 3 

(0.0152x + 0.468y)(3.020x - y) 4 

11.91 

6 

(. x 2 — 10 y 2 )x 3 

(5x + y)y 4 

11 

8.12 

(1.473x 2 + 5.397 xy + 2.867 y 2 )(x - 0.6732y) 3 

(-2.301x + 1.875y)(0.6732x + y) 4 

8.88 


The upper left quintic / in A( 3i1j1 ) is closest to h, at distance 3.0215666805997121633 1 / 2 , with 
triple root (—0.78519451639408253233 : 1). The lower right quintic g in A( 4j1 ) is closest to h, 
at distance 8.8808277614588859783 1 / 2 , with quadruple root (-1 : 0.67321557299682647408). 
We could have guessed the decomposition h — f + g in the fourth row from the input (J32J). 
This one is indeed a critical point, but it is neither primal optimal nor dual optimal. 

Our five solutions to (l3Tf were found by using the matrix that was introduced in Section 5: 


( 1 

5 

0 

-10 

0 

1 \ 

—s 3 

3 s 2 t 

—3 st 2 

t 3 

0 

0 

0 

—s 3 

3 s 2 t 

—3 st 2 

t 3 

0 

0 

0 

-s 3 

3 s 2 t 

—3 st 2 

t 3 

f 4 

4st 3 

6s 2 t 2 

4s 3 f 

s 4 

0 

V 0 

f 4 

4sf 3 

6 s 2 t 2 

4s 3 f 

s4 / 


The triple (resp. quadruple) roots (s : t) E P 1 of / (resp. g) are the roots of the binary form 
det (M 5 , 3 (s, t)) = (s 2 + t 2 ) 6 • (s 5 — 10s 3 t 2 — s 2 t 3 T 5st 4 ) = — (s 2 + t 2 ) 6 • (L ( 2 ^(/i)) (s, t). 


We compute the five real roots numerically, and at each of them we compute the left kernel 
of M. 5j3 (s,t). The result of that computation is precisely the list of five pairs (/,<?) above. 

This method scales well for hooks A = {l n-a , a}, ffere, the special ED degree is always n, 
and Lemma [5.21 furnishes the minimal polynomial for the desired a-fold root of f E A\. The 
matrix A i n , a (s,t) represents our task as a homogeneous polynomial eigenvalue problem, and 
this makes it amenable to well-developed methods of numerical linear algebra; see e.g. [7j. 
For an illustration we fix n = 15 and a = 6. The data point is the binary form 

h(x,y) = X^=o ^') u i xl y lb ~ l i with randomly chosen coefficient vector 

(u 0 ,..., u 15 ) = (20, -17, 3,16,12,14, -16, -5, 7, 8, -13, 5, -13, -16, 7, -11). 
Among forms / of degree 15 with a root of multiplicity 6, the following is the closest to h: 

Ei=o Ci) v i x i V 15 ~ i = 10“ 6 (® - 8.70886y) 6 (131903x 9 +11375.9x 8 y- 552.90hq/ 8 +45.8419y 9 ), 

{v 0 , vi , iq 5 ) = (20, -17, 3,16,12,14,-16, -5, 7.04, 8.19, -12.17, 8.09, -3.78,1.42, -0.46,0.13). 
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All 15 critical points for this optimization problem are real because L^\h) is real-rooted for 
our h. Its 15 roots have the form (s : 1) where s£i They appear in the first column of 


root s of L^ 5 \h) 

distance 2 

local 

8.70886 

86791 

min 

3.70567 

111796 

min 

2.19850 

163470 

max 

0.05736 

476068 

? 

-0.38870 

550056 

? 

-3.49092 

564363 

min 

-5.71229 

565936 

min 

-0.22118 

657621 

? 

1.25359 

723240 

max 

0.25811 

727831 

? 

0.48187 

774941 

? 

0.80694 

934884 

max 

-0.68808 

1058800 

max 

-1.67383 

1150260 

max 

-1.06515 

1256200 

max 


By computing the left kernel of A4 ^(s, t) at each root, we find the 15 decompositions (IHTTh 
The squared norms ||g || 2 of the dual solutions g = h — f are listed in the second column. 
So, the first row gives the optimal solution for A = {l 9 , 6 }, while the last row gives the 
optimal solution for the dual problem g = {l 4 ,11}. Local optima are indicated in the third 
column. There are four local optima on A^ (marked with “min”) and six local optima on A^ 
(marked with “max”). These were certified using the signature of the Hessian of the distance 
function. The local nature of the other five points cannot be decided with the second-order 
criterion because their Hessians are singular, in both the primal and the dual formulation. 

Our optimization problem is more challenging when A is not a hook. A small interesting 
case is the partition A = (3, 2 ). Let us see what happens here if we take the same input h as 
in ([32]) . We seek a binary form / E A( 3i2 ) with a triple root and a double root that is closest 
to h. The desired decomposition (I3T1) on the conormal variety Con( 3)2 ) now takes the form 

h(x,y ) = a(tx — sy) 3 [yx — uy) 2 + ((3x + r yy)(sx + ty) 4 + 5{ux + vy) 5 . 

This means that the vector (—1, a, (3, 7 , 5) lies in the left kernel of the 5 x 6 -matrix 

/ 1 5 0 -10 0 1 \ 

— s 3 u 2 2s 3 uv+3s 2 tu 2 —6s 2 tuv—s 3 v 2 —3st 2 u 2 3s 2 tv 2 +6st 2 uv+t 3 u 2 —3st 2 v 2 —2t 3 uv t 3 v 2 
0 t 4 4 st 3 6 s 2 t 2 4s 3 £ s 4 

t 4 4s£ 3 6 s 2 t 2 4s 3 £ s 4 0 

\ v 5 5 uv 4 lOtt 2 ?; 3 10 u 3 v 2 5 u 4 v u 5 ) 

All 5x5 minors of this matrix must be 0. However, the ideal of 5 x 5 minors has some 
extraneous associated primes that must be removed. This is analogous to the factor s 2 + 1 2 
in the hook case, but more complicated, so we do not pursue that primary decomposition. 
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Instead we simply work directly with the squared Euclidean distance function 

D =|| h — a(tx — syf(vx — uyf || 2 , 


and we solve the following system of polynomial equations in five unknowns: 

dD dD dD dD dD 

= = =-^— = o and sv-ut t o. 

da os dt du dv 

This has 20 complex solutions, while the ED degree of A( 3j2 ) is 21. One checks that h lies in 
the ED discriminant [SJ §7]. Only 4 of the 20 solutions are real. Setting t — v — 1, they are 


a 

s 

u 

distance 2 

1.817238 

-0.266252 

-0.265424 

1.815280 

0.673272 

-3.020572 

3.131909 

-0.785143 

-1.316853 

0.274673 

-0.387044 

1.428712 

7.724678 

8.643701 

12.017703 

13.105926 


The worst critical point is given in the bottom row. The corresponding quintic is 


/ = 1.81528a; 5 - 0.911262 x A y - 5.15521 x 3 y 2 + 0.0137459a;V + 4.34203 xy A + 1.79341a/ 5 . 


ED duality gives us the optimal solution g — h — f for the dual problem on (A( 3]2 )) v , namely 

g = -0.81528a; 5 + 0.911262a; 4 ?/ - 4.84479a;V - 0.0137459a;V + 0.657973a;?/ 4 - 0.793414?/ 5 

{fix + 7 |/) (—0.785143a; + yf + 5(1.428712a; + yf 

for suitable real constants /3, 7 , and S, which can be found by solving a linear system. 

The form g is very interesting for the application described in Section 4. Recall that the 
generic tensor rank for binary quintics is 3. The unique rank 3 decomposition of the given 
quintic h is complex. It is shown on the left in f[32|l . One might therefore ask for a best 
approximation to h that has real rank 3. This question is not well-posed because the set 
77-5 from Section 4 is not closed. Instead one should ask for the closest binary quintic in the 
closure of 1Z 5 , i.e. among those whose border rank equals 3. That optimal quintic must be our 
g, because it is a critical point on real rank boundary d^giTZ^), which equals the little apple 
hypersurface (A( 3 2 )) v . The fact that g has border rank 3 is verified by the representation 

g = lim ^((-0.785143 + e@)x + (1 + eq)yf - ^(-0.785143a; + yf + <J(1.428712x + y)^j . 

This computation offers a concrete illustration of Theorem 14.11 and Example 14.21 


Acknowledgements. 

We thank Michael Burr for a conversation in March 2015 that ignited our interest in the ED 
problem for multiple root loci. We are grateful to Greg Blekherman for his help in getting 
Theorem 14.II into final shape. This project was carried out in the summer of 2015, when both 
authors visited the National Institute of Mathematical Sciences, Daejeon, Korea. Hwangrae 
Lee was supported by the National Research Foundation of Korea (MSIP 2011-0030044). 
Bernd Sturmfels was supported by the US National Science Foundation (DMS-1419018). 


22 












References 


[1] A. Abdesselam and J. Chipalkatti: Brill-Gordan loci, transvectants and an analogue of the 
Foulkes conjecture , Adv. Math. 208 (2007) 491-520. 

[2] A. Abdesselam and J. Chipalkatti: The bipartite Brill-Gordan locus and angular momentum , 
Transform. Groups 11 (2006) 341-370. 

[3] G. Blekherman: Typical real ranks of binary forms, Found. Cornput. Math. 15 (2015) 793-798. 

[4] G. Blekherman, P. Parrilo and R. Thomas: Semidefinite Optimization and Convex Algebraic 
Geometry, MOS-SIAM Series on Optimization, 13, 2012. 

[5] J. Chipalkatti: On equations defining coincident root loci, J. Algebra 267 (2003) 246-271. 

[6] P. Comon and G. Ottaviani: On the typical rank of real binary forms, Linear and Multilinear 
Algebra 60 (2012) 657-667. 

[7] J-P. Dedieu and F. Tisseur: Perturbation theory for homogeneous polynomial eigenvalue prob¬ 
lems, Linear Algebra Appl. 358 (2003) 71-94. 

[8] J. Draisma, E. Horobel;, G. Ottaviani, B. Sturmfels, and R. Thomas: The Euclidean distance 
degree of an algebraic variety, Found. Comput. Math., to appear. 

[9] L.M. Feher, A. Nemethi and R. Rimanyi: Coincident root loci of binary forms, Michigan 
J. Math. 54 (2006) 375-392. 

[10] I.M. Gel’fand, M.M. Kapranov and A.V. Zelevinsky: Discriminants, Resultants and Multidi¬ 
mensional Determinants, Birkhauser, Boston, 1994. 

[11] D. Hilbert: Singularitaten der Diskriminantenflache, Math. Annalen 30 (1887) 437-441. 

[12] A. Iarrobino and V. Kanev: Power Sums, Gorenstein Algebras, and Determinantal Loci, Lec¬ 
ture Notes in Mathematics, 1721, Springer-Verlag, Berlin, 1999. 

[13] E. Kaltofen, J. May, Z. Yang, and L. Zhi: Approximate factorization of multivariate polyno¬ 
mials using singular value decomposition, J. Symbolic Comput. 43 (2008) 359-376. 

[14] G. Katz: Plow tangents solve algebraic equations, or a remarkable geometry of discriminant 
varieties, Expo. Math. 21 (2003) 219-261. 

[15] S. Kurmann: Some remarks on equations defining coincident root loci, J. Algebra 352 (2012) 
223-231. 

[16] J. McNamee: Numerical Mathods for Roots of Polynomials, Part I, Studies in Computational 
Mathematics, 14. Elsevier B.V., Amsterdam, 2007. 

[17] E. Miller and B. Sturmfels: Combinatorial Commutative Algebra, Graduate Texts in Mathe¬ 
matics, Springer Verlag, New York, 2004. 

[18] F. Napolitano: Topology of complements of strata of the discriminant of polynomials, Comptes 
Rendue Acad. Sci. Paris, Serie I, 327 (1998) 665-670. 

[19] L. Oeding: Hyperdeterminants of polynomials. Adv. Math. 231 (2012) 1308-1326. 

[20] G. Ottaviani, P.-J. Spaenlehauer and B. Sturmfels: Exact solutions in structured low-rank 
approximation, SIAM J. Matrix Anal. Appl. 35 (2014) 1521-1542. 

[21] B. Sturmfels: The Hurwitz form of a projective variety, arXiv: 1410.6703. 

[22] J. Weyman: The equations of strata for binary forms, J. Algebra 122 (1989) 244-249. 

[23] Z. Zeng: Computing multiple roots of inexact polynomials, Math. Comp. 74 (2005) 869-903. 


23 


[24] Z. Zeng: Regularization and matrix computation in numerical polynomial algebra , in Approxi¬ 
mate Commutative Algebra, 125-162, Texts Monogr. Symbol. Comp., Springer, Vienna, 2009. 


Authors’ addresses: 

Hwangrae Lee, Pohang University of Science and Technology, Korea, meso@postech.ac.kr 
Bernd Sturmfels, University of California, Berkeley, USA, bernd@berkeley. edu 


24 


